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Abstract 

In this work we study a modified theory of gravity that contains up to fourth order spatial deriva- 
tives as a model for the Hofava-Lifshitz gravity. The propagator is evaluated and, as a result, it is 
obtained one extra pole corresponding to a spin two nonrelativistic massless particle, an extra term 
which jeopardizes renormalizability, besides the unexpected general relativity unmodified propaga- 
tor. Then, unitarity is proved at the tree-level, where the general relativity pole has shown to have 
no dynamics, remaining only the two degrees of freedom of the new pole. Next, the nonrelativistic 
effective potential is determined from a scattering process of two identical massive gravitationally 
interacting bosons. In this limit, Newton's potential is obtained, together with a Darwin-like term 
that comes from the extra non-pole term in the propagator. Regarding renormalizability, this extra 
term may be harmful, by power counting, but it can be eliminated by adjusting the free parameters 
of the model. This adjustment is in accord with the detailed balance condition suggested in the lit- 
erature and shows that the way in which extra spatial derivative terms are added is of fundamental 
importance. 
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I. INTRODUCTION 



Einstein's General Relativity (GR) generalizations due to the addition of extra derivative 
terms was proposed a long time ago jl|, I2I. Such modifications became interesting in the 
context of quantum gravity, where Einstein's theory is known to be perturbatively nonrenor- 
malizable by power counting jsl. It has been verified that, in the presence of extra terms 
containing products of the curvature tensor, the theory turned out to be renormalizable. 
However, such modification introduced pathologies into the theory which loses unitarity j4|. 

More recently, in an attempt of constructing a perturbatively renormalizable and unitary 
gravitational theory, Hofava jsl proposed modifications of GR via extra terms with only 
spatial derivatives, introduced in a chosen spacetime foliation. The foremost argument for 
such proposal lies in the fact that, the gravity propagator of the linearized theory would 
behave like 

(1.1) 



w2 _ p _ a2(A;2)2 a,{k^y 

a2, ■ ■ ■ , ttz being coupling constants, k'^ = {u,k) the four-momentum of the graviton and 
z > 1 a parameter associated with the highest order of spatial derivatives. The absence of 
higher order time derivatives may transform the propagator as the one shown in (11. ip . with 
only simple poles in cj^. It has been verified that the theory proposed in ^ is not so simple 
and that new degrees of freedom, among other illness, are present jo - 14|. 

In this work we will be interested to study the exact form for (II. ip . To this end, we 
propose a prototype for the theory in [Si]. The model in question, restricting ourselves to 
z = 2, is described by the action 

S = l f dt [ (fx^g + ^ f^^^JJ^ ^ (1.2) 

defined in the foliation = 3ft x a. There, spacetime indices fi, u, ■ ■ ■ run from to 3, 
while i, j = 1,2,3 are the indices that label spacial coordinate on a. Also, = 16ttG, 
G is the Newton's constant, and, all over the paper, c = h = 1. The spacetime signature 

we are dealing with is ( h ++), R^^^'^p iR-^%i) ^^e Riemann curvature tensor in 4 (3) 

dimensions for the metric g^^^, (gij). It will be convenient, in this spacetime signature, to 
define = d^vi'}'' - 9^rg" + ■ ■ ■ , while i^^fj = i?^^. Similar definitions are valid for 

the 3-dimensional curvature tensor. 
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Notice that, to be as general as possible, we could have incorporated the term R^'^^^^''''R[^Ij 
into (II. 2p . However, in the 3- dimensional case the Weyl tensor jl| is identically zero. As a 
consequence, the aforementioned term R^^^^^'^^R^'^l^ can be totaly cast as a combination of 



/j(3)2 g^j^^ R^'^^'^^ Rf^ (see, for instance, the appendix in 15|). 

In the next section we will dedicate ourselves to the computation of the propagator of 
f ll.2p in the weak field approximation. Section IIIII contains a systematic study of unitarity 
at the tree-level. At this level, we show that the dynamic of the theory is governed by a 
pole that corresponds to a nonrelativistic spin two massless particle. The obtention of the 
semiclassical nonrelativistic potential for a boson-boson scattering process, via gravitational 
interaction, is done in Section llVi Section |V] contains the conclusions. 

II. THE PROPAGATOR 

Let us perform the linearization of (11. 2p in the weak field approximation 

g^lu'^r]^^, + Kh^^, (2.1) 

whose inverse is g^^ ^ rj^^ — K,h^^ . The background metric is rj^y = diag{ — l, 1, 1, 1) and the 
hfj_u are the gravitational field fluctuations. By collecting terms up to second order in h^^, 
(USD gives 



S=- I d^x -h^^d^h"'' --hd^h + r + 2aR^^^^ + 2(3R^^^'^R!:f 
— 2 J 2 4 J 



(2.2) 



where d"^ = d^d^, simple underline means first order in h, double underline means second 
order in h and so on. The trace is h = rj^'^h^j^ while 



— 'I L-al3 

= do^W - ^d^h . (2.3) 

Clearly, the action (12. 2p is invariant under the gauge transformations Sh^^i, = —d^C,u " 
dy^^ = —2d(^fj_C,u) for any arbitrary C,ii- Then, in order to evaluate the propagator, we may 
choose the de Bonder gauge, = 0, by introducing into the action the following term: 

S,f = ~ jd'xT^T^. (2.4) 
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Now, we may rewrite (12. 2p as 



(2.5) 



where the operator O possesses the symmetries Of^u,ai3 = Oai3,fiu = Oi,^^^ap- It is convenient 
to separate O into 



where in momentum space, with = {oj, k), 

k"^ 9 - \ 



k'^V^iuVaf} + (1 - >^)Vmvkp)k^)) 



1-A 



ij^fipkakp + Tjapk^kp 



(2.6) 



(2.7) 



whereas 0^^^ is equal to zero whenever one of its indices is time-hke and 

^liki = Tik'^^tjM + (2a + (3)kikjkikk + {-2a - ^)fc^ {Sijkkh + Skikikj) 
■"I I 

+(2a + -)k'^5ijSki - pP5^^i(^kh)kj)) . 



(2.8) 



In the above expressions we have set fc^ = — w^ + fc^, S^y^^p = VtJ'{ciVi3)iy used the convention 

A{{fi{aBi3)u)) = {I / 2) {A ^(^a 3/3)^ + y4y(Q,i?/3)^) . 

The difficulty of inverting ( 12. 6p comes from the noncovariant form of the operator O. Our 
plan is to ffist separate the components with pure spatial indices from the rest. This may 
be achieved by defining 



A 



ij,kl 



^ij,kl + ^ij,kl ■ 



(2.9a) 
(2.9b) 



Componentwise, 



^00,00 



Aoi,kl 



M.J 



4 ' 



2 ^ 
1-A 



A 



X 



_ p _ (1 _ X)co' 



1-A 

~l~ „ kiki 



Oiflj 

^5ij [e + (1 - \)uj'] - , 



(2.10a) 
(2.10b) 
(2.10c) 

(2.10d) 
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and zero otherwise. We now may write 



(2.11) 



From the inverse equation 



one obtain 



o 


-100,00 


o- 


1 00,m?i 


o 


-lOi.OO 


&- 


1 Oi,mn 




1 kl,mn 



— (i-Aoo.,a-^^^'°°) 

in nn V / 



^00,00 

A. 



00,kl Q-lkUmn 



^00,00 
1 



1 



-lfc«,00 



ij 



By the definition of Mjj in (I2.10dl) . it is straightforward to get 
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F + (1 - A)w2 
In the last fine of f l2.13p we defined 

^00,A:/^ij,00 



a 



ij,kl 



^i-i.ki ( + 



'ij,kl 

kikjkikk 



k^ 

Sijk^ki + dkikikj 



^00,00 

2 2 

(2a + /3)P + 



:4 , (1 - >^?k^ 



(1 - Xfu^k^ 



AP AP[A;2 + (1-A)a;2] 



A;2 

(1 - A)2a;2g2- 

AF 
(1 - A)2a;2^2- 

AA;2 



(2.12) 

(2.13a) 
(2.13b) 

(2.13c) 

(2.13d) 
(2.13e) 



(2.14) 



45((i(fcA;/)%)) 



2- A 



1 /2 - A 



(2" + t)^ ■ AFV 2 



+ 



A;2 



4 



4 



(2.15) 



4[A;2 + (1 - A)a;2] 

The appendix [X] contains all the necessary tools to invert Cij^ki- The calculation is tedious 
but straightforward. We use the Barnes- Rivers operators given in (lAip . together with (lA4p . 
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5]), and, as a last step in the equation so obtained, the equahties in (lASP and (1A9P are 



employed. We then quote the result 

P 



(2P^ + 4P°).. 



ij,kl 



ij,kl 



k' 



^i^-(8a + 3/3)^P° 



(2pi + 4po) 



ij,kl 



ij,kl 



(2.16) 



So far we have obtained the part of the propagator with pure spatial indices. Equations 
f l2.13ap . ( ]2.13bp . (I2.13cp and (12.13dp enable us to get the remaining terms, so that 



2pi + 4-po _ -po _ 3p.o _ p.0 
+ 



-XP 



'^^%^al3'^Pij,kl 



-P 



-P + jSk^ 
8a + 3/3) Q^^^ap , 



(2.17) 



where 



Q 



Qij,kl — 'Pij,kl ; 
^00,00 = ^ , 
0" 



,mn 



k* 



Q 



mn,00 



0, otherwise. 



(2.18) 



The propagator in f l2.17p has two poles besides the term Q. The first line of that equation 
is just what one gets from pure GR, and corresponds to the massless pole —k'^ = u'^ — k'^ = 0. 
By looking at the second line in (I2.17p . we notice that this pole, for nonvanishing /3, gains a 
correction that is proportional to in its spatial indices sector. In the next section we will 
analyze the contribution of this correction to the dynamic of the theory, at the tree-level. 
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Besides this just mentioned pole, there is a new pole corresponding to a massless spin 2 
particle which obeys a nonrelativistic dispersion relation 



(2.19) 



For this pole to have physical meaning either /3 < 0, or, otherwise, there will be a limit in 
the particle momentum (/3 > ^ /c^ < 1//3). This is the expected pole in the propagator 
we wrote in (II. ip that may improve renormalizability. Nevertheless, the last term in ( 12.171) . 
proportional to Q, clearly spoils renormalizability, unless we set 8a + 3/3 = 0. Such choice 
of the parameters a and f3 is in accordance with the detailed balance condition introduced 



by Hofava in 5j. This can be seen by taking a = — 1/3 so that the extra spacial derivative 



terms in the action (11.21) furnishes 



13 / rf^x. 



/3 / d'^Xy 



-g (^J_Rm + i^f^ 

-5W[q] 5W[q] 

'9—r yij,kl~i ) 

dqij oqu 



(2.20) 



where 



QijM — (li{kqi)j — 2^ij(lkl 



is the inverse of the Weyl metric = g^C^gOi _ 



kl 



qij and the 3-dimensional action W is given by 



W[q] 



d^Xy/qR^^^ . 



(2.21) 

qij = Qij while, g*-^ is the inverse of 

(2.22) 



It is worth mentioning that the detailed balance condition plays an essencial role regarding 



renormalizability in the full Hofava theory, as shown in 



16] . In the present case it justifies 



the removal of a bad behaved term which would spoil the renormalization of the model. 



III. TREE-LEVEL UNITARITY 



We can readily check from (I2.17P that the linearized theory described by (11.21) does not 
contain tachyons. This is already an improvement compared to the modified theories of 



gravity with higher time derivatives jj, ll7| that, beyond having tachyons, are not unitary 
even at the tree-level. 
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In this section we are going to examine unitarity for the model (11. 2p at the tree-level. 



To this end, it is enough to study the residue o: 
(I2.17P with an arbitrary conserved current 17|, 



the poles when we saturate the propagator 



T^^^ = akf'k" + 6Fr + c^ye^e^ + 2dkH''^ + 2eJ^^k''^ + 2f^ei^k''^ . (3.1) 

Here, T^^^ has been arbitrarily expanded in terms of the linearly independent four-vectors 
k^ = {00, k), k'^ = {—u,k) and the orthonormal graviton polarization vectors = (0,e,;), 
x,y = 1,2 {e^ ■ k = 0) with the coefficients a, b, c^y = Cyx, d, and f^. Conservation implies 
T^-'k^ = 0, so that 

ak'^ + d{uj^ + P) = , (3.2a) 
dk^ + b{u^ + p) = , (3.2b) 
e,P + /,(u;2^P) =0, (3.2c) 
ak* + b{uj^ + P)2 + 2dk\uj^ + P) = . (3.2d) 

Equations (13. 2p and (13. ip enable us to write the amplitude 

r^ll _ ^22n2 , 4/ (k^ - 3fc2p + 

T^'O-'sT'^^ = ^ ^ ^ - 16(8a + 3/3)^^0^^ , (3.3) 

-k' + Pk^ [u^ + k^y 

where we used the identities TVT = ioT V = and 

= {c''-c^^Y + A{c'y>0. (3.4) 

The last term in (13. 3p does not correspond to a pole, so, there is no particle associated 
to it. But, it clearly prejudices the divergent behavior in perturbation theory unless we fix 
accordingly the constants a and /3, i.e., 8a + 3/3 = 0. In this situation, we conclude that, 
at least at the tree-level, the theory is unitary. The cancellation of the term corresponding 
to the pole at —k"^ = in (13. 3p implies that this particle has no dynamics. Yet, the pole 
at —k"^ + iSk"^ = with positive residue is a physical particle with two degrees of freedom 
{Pfjki has only two independent indices). In other words, apparently, the disappearance of 
the pole — A;^ = in favor of the pole — A;^ + (ik"^ = means that the graviton has become 
a nonrelativistic particle represented by this pole. Notice that the residue of this modified 
pole is the same as the graviton pole when a = (3 = 0. 
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FIG. 1: Second order Feynman diagram of two bosons interacting via gravitational exchange. 

IV. SEMICLASSICAL POTENTIAL IN THE NONRELATIVISTIC LIMIT 

Once the propagator of the theory is obtained, we can analyze the scattering process 
of two particles interacting gravitationally within this modified theory. This enable us to 
evaluate the effective low energy potential due to the gravitational interaction of two identical 
massive bosons particles of zero spin described by the Lagrangian density 

L = [—d,_iipd^(p* — m^ifip*^ 

~ ^matter + ; (4-1) 

up to first order in h, where 

^matter = -Of^ipd^ip* - Ul^ipip* , (4.2a) 

The amplitude for this scattering process, as illustrated in Fig. [H is given by 

where k = p' — p = q — q' is the momentum transfer and the vertex Vf^i^{p,p'), drawn in 
Fig. [21 is 

Vt^iyiP^P) = -2 [^Pif^Piy) - Vt^^ [P-P - m^)] , (4.4) 

with p-p' = p^p'^. For simplicity, we will restrict our calculation to the center of momentum 
reference frame (p^ = {E,p), p'^ = {E,p'), q^ = {E,—p), and g'^ = {E,—p')), such that 
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FIG. 2: Boson-boson interaction vertex. 



the amplitude (14. 3 P can be cast as 



K 



T 



2F 



[{q ■ p + q ■ p') {q ■ P + q' ■ p') — {q ■ q' + {p ■ p' + m^) — 2m^] 



+ - 



{q-p + q-p'){q'-p + q"-p')--q-{q + q ')p ■ {p + P ') 



fc2 A;2 - (3k' 



-(8a + 3/3)^ {2pY' + p-p' + m^Y 



2A;2 



[{q ■ p + q ■ p) {q ■ p + q ■ p) ~ {q ' q + t^'^) {p ■ p + — 2m^] 



H — r(P'5'+"^) ^ — 

4 ^ ' VP P 



(4.5) 



The effective potencial for the just calculated scattering amplitude is obtained from the 
Fourier transform 



V{x) 



1 



8m2 } (27r)3 

by inserting the low energy limit of (14. 5p . that is to say. 



TNRexp{—ik ■ x) , 



Tnr 



2 4 

Km 



A;2 



(8a + 3/3)fi:W. 



Collecting all the above results one gets 

V{x) = 



\x\ 



-2{8a + 3f3)Gm^6^^\x) 



(4.6) 



(4.7) 



(4.8) 



As one can observe, the nonrelativistic limit (14. 7p eliminates the pole — A;^ + /3A;^ = 0. The 
pure GR sector obviously reproduces Newton's potential. But this is not all. The extra non- 
pole term in (I2.17P is the responsible for the appearance of a Dirac's delta in the potential. 
This is similar to what happens in QED in the calculation of the nonrelativistic effective 
potential of two electrons exchanging a photon. In that case, this extra delta term, also 
known as Darwin's term, is interpreted as quantum fluctuations in the electron's position 
due to its position indeterminancy. 
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V. SUMMARY AND CONCLUSIONS 



In this paper we have studied the propagator of a Horava-Lifshitz hke theory with quartic 
spatial derivative terms. This propagator has two poles, one corresponding to the GR gravi- 
ton pole and other also corresponding to a spin 2 massless but of nonrelativistic character. 
Besides the improvement brought by the spatial quartic terms there are some points which 
deserve mentioning. Firstly, the presence of the unmodified GR sector is potentially dan- 
gerous. Nevertheless, by coupling the model to a scalar field through a conserved current, 
we verified that at tree-level the residue of the GR pole is zero and that the only excitation 
corresponds to a spin 2 nonrelativistic particle. In that situation, we also showed that the 
theory is free of ghosts and tachyons. Even though at the tree-level the unmodified GR 
pole has no dynamics, we cannot assert that it will not contribute to higher order virtual 
processes. Moreover, the propagator also possess a non-pole term that by power counting 
clearly prejudices renormalizability . In fact, such term increases the ultraviolet divergence 
and may spoil the theory as a whole. The solution for this problem passes by a choice of 
the arbitrarily inserted constants a and /3 and, as showed, is in the class of extensions which 
satisfies the detailed balance condition. 

The effective low energy potential for scattering process involving two massive bosons 
that interact via this higher spatial derivative theory was also computed and, as a result, we 
obtained Newton's potential plus a Dirac delta in position, i.e., a Darwin-like term. This 
term is well known in QED from the evaluation of the nonrelativistic potential obtained 
from electron-electron scattering and is related to the indeterminancy in electrons' position. 
In the quantum gravity case, this imprecision in the boson position is not present in pure 
linearized GR and have appeared in the present case labeled by the constants a and /3. 

The relevant information we have obtained with this model is that higher order spatial 
derivatives do not ensure that the GR propagator will be modified as expected, with higher 
order in momentum fc^. In fact, within the modifications worked in this paper, we have 
reobtained the unwished GR graviton pole plus the desired term like (II. ip which may improve 
the ultraviolet behavior of the quantum theory. We also showed that a bad ultraviolet 
behaving term has appeared and that only within a specific combination of the constants a 
and /3 it can be eliminated. This, in fact, shows that the addition of higher order spatial 
derivatives is not enough to warrant renormalizability and that the way one introduces such 
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extra terms is crucial when renormalizability is at stake. 
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Appendix A: Barnes-Rivers operators 



The 3-dimensional symmetric Barnes-Rivers operators 



17 



19] are given by 



Pij,ki = 26'((i(fcW,)j)) , (Ala) 
Pij,ki = (^i{kOi)j - -^dijdki , (Alb) 



P?iM = , (Ale) 



Pij,ki = ^ij^ku (Aid) 
Pfj ki = Gij^ki + ^ijdki , (Ale) 



where the projection tensors 



Oij - Sij ^ , Uij - -jr^ (A2) 

have been defined. Such operators obey (using AB in the place of A^^'^^Bki mn to the con- 
tractions) Pipi = P\ P^P^ = P\ P'^P'^ = po, popo = po, popo = {D- l)(po + PO), 
popo ^ popo ^ peu^ popo ^ popo ^ pu^e^ together with Pj^^^^^i = Oijooki and P^'j^ = UijOki- 
Any other contraction is found to be zero. Those operators also obey the identities 



5,,4i = (2P° + P° + P%,H, 

f^%(kki)kj)) = i2P' + AP%j,ki , 
^ {5i^kkh + 5kihki) = (P° + 2P°) 



-^{kikjkkh) = Pij ki ■ 
k* 



ij,kl 



(A3a) 
(A3b) 
(A3c) 

(A3d) 

(A3e) 
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The above tools enable one to write an arbitrary symmetric operator O as 



(A4) 



whose inverse, if exists, will be 



1 



XqP^ + XoP^ - XqP 



Xi X2 XqXo — 2Xq 

In 4-dimensions, the Barnes-Rivers operators may be written as 



0). 



(A5) 



' iJ.u,a.p 







(A6a) 

(A6b) 

(A6c) 
(A6d) 
(A6e) 



Now, the projection operators are defined by 

kfjku 



n 



(A7) 



' k'^ ' 

Eventually, it will be convenient to relate the Barnes-Rivers operators in three and four 
dimensions. When only the spatial indices are being treated, it is possible to obtain the 
identities between the P's and T^'s as 

— ^ - -v^ + — p° - — : 

fc2 2 2A;4 2A;2 



p2 



ij,kl 



ij,kl 



pO 

ij,kl 



pO 

ij,kl 



pO 
ij,kl 



k'^ k^ / ij,kl 

2 2k^ 2P 



ij,kl 



k' 



k^ 



' ij,kl 1 



^2 ^4 



(A8a) 

(A8b) 

(A8c) 

(A8d) 
(A8e) 



Such equations enable one to get 
^2p2 _ _ _ p^ 



ij,kl 



A;2 



A;2 



(A9) 



ij,kl 
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